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ABSTRACT : 

Experimental  results  for  the  convective  heat  transfer  and  flow 
friction  characteristics  of  plate-fin  type  heat  transfer  surfaces 
are  presented  for  five  surfaces.   Four  surfaces  had  a  sinusoidal  fin 
geometry.   Two  of  these  were  fabricated  of  stainless  steel,  two  of 
nickel.   One  of  the  nickel  surfaces  had  perforated  nickel  fins.   The 
fifth  surface  was  fabricated  of  stainless  steel  and  had  a  triangular 
fin  geometry.   The  heat  transfer  data  were  obtained  by  the  single-blow 
transient  testing  technique  and  by  the  cyclic  testing  technique.   The 
two  techniques  compliment  each  other  in  extending  the  range  of  testing 
capability. 
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NOMENCLATURE 

English  Letter  Symbols 

A      Matrix  total  heat  transfer  area      sq  ft 

A„     Matrix  minimum  free  flow  area       sq  ft 

c  ;* 

Ar  Matrix  total  frontal  area  sq  ft 

As     Matrix  solid  cross  sectional  area    sq  ft 
available  for  thermal  conduction 

At      Matrix  solid  cross  sectional         sq  ft 
conduction  area  considering 
effect  of  matrix  material 
perforations 

ag     Fin  thickness  ft 

C      Jet  contraction-area  ratio  for       dimensionless 
circular  tube 

c      Specific  heat  (gas)  at  constant      Rtu/(lbm  deg  F) 
pressure 

cs     Matrix  material  specific  heat        Etu/(lbm  deg  F) 

Cf     Flow  stream  capacity  rate  (mcp)      Btu/(hr  deg  F) 

Cg     Matrix  heat  capacity  (Wscg)  Btu/(deg  F) 

D^     Hydraulic  diameter  of  internal       ft 
passage 

E      Friction  power  expended  per  unit     hp/sq  ft 
of  surface  heat  transfer  area 

G      Matrix  flow  stream  mass  velocity     lbm/(hr  sq  ft) 
(m/Ac) 


g      Proportionality  factor  in  Newton's    32.2  (lbm  ft)/(lbf  sec  ) 
Second  Law 

h      Thermal  convection  surface  heat      Btu/(hr  sq  ft  deg  F) 
transfer  coefficient;  heat  transfer 
power  per  unit  area  per  degree 
temperature  difference 

K      Loss  coefficient  for  sudden  flow     dimensionless 
c 

contraction  at  matrix  entrance 


K,     Momentum  velocity-distribution 
coefficient 

Ke     Loss  coefficient  for  sidden  flow 
expansion  at  matrix  exit 

k      Unit  thermal  conductivity 

ks     Matrix  material  thermal  con- 
ductivity 

ky     Thermal  conductivity  in  radial 
direction 


Total  matrix  flow  length 

Mass  flow  rate 

Pressure 

Matrix  porosity  (A  /A,  ) 

Heat  transfer  rate 

Universal  gas  constant  (53.35  for 
air) 

:„     Hydraulic  radius  (ACL/A^) 


L 
m 
P 
P 

q 

R 


dimensionless 

dimensionless 

Btu/(hr  sq  ft  deg  F/ft) 
Btu/(hr  sq  ft  deg  F/ft) 

Btu/(hr  sq  ft  deg  F/ft) 

ft 

lbm/hr 

lbf/sq  ft 

dimensionless 

Btu/hr 

(ft  lbf)/(lbm  deg  R) 

ft 


m 


V. 


W. 


Relative  shape  factor 

Temperature 

Flow  velocity 

Matrix  volume 

Matrix  material  volume  (corrected 
for  effects  of  perforations) 

Matrix  mass 

Distance  from  matrix  inlet  in 
flow  direction 


Subscripts 


atm    Local  atmosphere 
ave    Average 
f      Fluid  (gas,  air) 
i      Initial,  inlet,  individual 
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dimensionless 
deg  F 
ft/sec 
cu  ft 
cu  ft 

lbm 
ft 


m  Mean,  matrix 

0  At  orifice 

s  Solid  (matrix  material) ,  static 

STD  Standard  (temperature  and  pressure) 

x  Local  conditions 

1  At  inlet  (upstream  of  heating  elements) 

2  At  matrix  inlet  (downstream  of  heating) 

3  At  matrix  exit 


Greek  Letter  Symbols 

P      Area  density:  compactness  (ratio  of 
matrix  total  heat  transfer  area  to 
core  volume-A/V) 

8      Diameter  ratio,  do/d 

A      Difference  or  change  (time,  distance, 
temperature) 

0  Time 

U  Fluid  viscosity  (dynamic) 

v  Specific  volume 

p  Density 

o      Ratio  of  free  flow  to  frontal  area 
(Ac/Afr) 

w0     Angular  frequency  of  temperature 
oscillation 


sq  ft/cu  ft 
dimensionless 


Sec,  hr. 
lbm/hr  ft 
ft3/lbm 
lbm/cu  ft 

dimensionless 

rad/sec 


Dimensionless  Groupings 


Amplitude  attenuation  parameter 

Fanning  friction  factor  (Ratio  of  wall  shear  stress  to  the 
fluid  dynamic  head).   This  factor  is  plotted  versus  N^  to 
illustrate  the  flow  friction  characteristics  of  a  matrix 
surface . 


j      Colburn  Factor  (NstN    )  ,  a  generalized  heat  transfer 

grouping.   This  factor  is  plotted  versus  Nr  to  illustrate 
heat  transfer  characteristics  of  a  matrix  surface. 

M  Flow  parameter  (mcp/Wcs  w0  ) 

N  Nusselt  number  (hD^/k  ),  a  heat  transfer  modulus 

Np  Prandtl  number  (ucp/k  ) ,  a  fluid  properties  modulus 

N  Reynolds  number  (D^G/y  ) ,  a  flow  modulus 

Ngt  Stanton  number  (h/Gcp  ) ,  a  heat  transfer  modulus 

Ntu    Number  of  heat  transfer  units  (hA/mcp) ,  a  heat  transfer 
parameter 

R      Ratio  of  matrix  outlet  temperature  amplitude  to  matrix  inlet 
temperature  amplitude  during  cyclic  temperature  variation 

X      Longitudinal  heat  conduction  parameter  for  solid  matrix 
material  (ksAs/mLcp). 

£  Distance  parameter  (x/L) 

Z  Time  parameter  (hA/Ws  Cs)  used  in  Single-Blow  analysis  (17) 

t  Time  parameter  (u©/L  )  used  in  cyclic  technique 

id  Frequency  parameter  (u>  L/u  )  used  in  cyclic  technique 


INTRODUCTION 

This  document  reports  the  results  obtained  during  FY66  and  FY67 
of  a  continuing  program  for  the  evaluation  of  the  heat  transfer  and 
fluid  friction  characteristics  of  compact  heat  exchanger  surfaces  at 
the  Naval  Postgraduate  School,  Monterey.   Results  obtained  in  FY65 
have  been  previously  reported  [1]  ,  [2]  . 

DESCRIPTION  OF  SURFACES 

The  five  surfaces  tested,  manufactured  and  provided  by  Solar,  a 
division  of  International  Harvester,  San  Diego,  had  basically  sin- 
usoidal or  triangular  fins  in  a  plate-fin  geometry.   Table  I  is  a 
summary  of  the  matrix  geometries.   Only  one  matrix  tested,  designated 
Solar  #3,  was  brazed.   All  matrices  had  a  frontal  cross  section  of 
approximately  3  and  1/16  inches  square  and  a  length  of  approximately 
3.  inches  in  the  flow  direction.   The  perforated  nickel  plate  fins  of 
the  Solar  #2  matrix  is  an  electro-deposited  metallic  sheet  of  integral 
structure  manufactured  by  Perforated  Products,  Incorporated,  and  is 
described  in  Appendix  D. 


1.  Numbers  in  brackets  refer  to  numbered  items  listed  in  References 


EXPERIMENTAL  METHODS 

To  obtain  the  heat  transfer  data,  two  testing  techniques  were 
employed:   the  single-blow,  transient  testing  technique  and  the  cyclic 
testing  technique. 

The  single-blow,  transient  testing  technique  is  completely  des- 
cribed in  reference  [2] ,  however,  briefly  this  method  consists  of 
heating  the  test  matrix  by  heated  air  to  a  uniform  temperature  (approx- 
imately 20. F  above  ambient  temperature  in  these  tests)  and  then  sub- 
jecting the  matrix  to  a  step  change  in  the  air  flow  temperature  to 
ambient  temperature.   The  air  temperature  downstream  of  the  matrix  is 
monitored  and  recorded  versus  time.   The  test  rig  used  in  these  tests 
is  shown  in  Figures  1-5,  and  described  in  detail  in  Appendix  C.   A 
series  of  fine  nichrome  wire  heaters  were  installed  upstream  of  the 
matrix  to  heat  the  air.   Turning  off  the  heaters  provides  the  step 
change  in  the  temperature  of  the  air  without  disturbing  the  flow.   By 
referencing  the  thermocouples  downstream  of  the  matrix  to  the  thermo- 
couples upstream  of  the  heaters,  the  initial  temperature  difference  can 
be  very  closely  controlled.   This  temperature  difference  is  continuously 
recorded  by  a  strip  chart  recorder  when  a  run  is  made.   This  recorded 
trace  has  the  distinct  advantage  that  there  is  no  transposition  of 
data  required  which  would  produce  increased  uncertainties.   Reduction 
of  the  data  follows  the  method  of  Locke  [6]  and  Howard  [4] ,  whereby  the 
Ntu  of  the  surface  can  be  evaluated  by  determining  the  maximum  slope 
of  the  fluid  temperature  difference  versus  time  curve  during  the  cooling 
transient. 


10 


The  cyclic  testing  technique  is  completely  described  in  reference 
[51.   However,  briefly  this  method  consists  of  producing  a  sinusoidal 
time  varying  temperature  of  the  steady  air  flow  to  the  matrix.   This 
air  inlet  tenmerature  and  the  time  varving  temperature  of  the  outlet 
air  downstream  of  the  matrix  are  monitored  and  simultaneously  recorded. 
The  test  rig  used  was  the  same  as  that  used  for  the  above  single-blow 
transient  tests.   A  low  frequency  sinusoidal  voltage  was  superimposed 
on  an  appropriate  d.c.  level  and  fed  to  the  nichrome  heaters.   Reduction 
of  the  data  follows  the  method  of  Bell  and  Katz  [5],  whereby  the  Ntu 
of  the  surface  can  be  evaluated  by  determining  the  ratio  of  output  to 
input  amplitudes  of  the  temperature  variations.   See  Appendix  A.   Again, 
recording  these  temperatures  direct lv  and  simultaneously,  requiring  no 
transposition  of  data,  reduces  the  data  reduction  uncertainties. 
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PRESENTATION  OF  RESULTS 

Heat  transfer  and  flow  friction  data  for  each  of  the  matrices 
tested  are  given  in  both  tabular  and  graphical  form,  TABLES  I  through 
IV  and  Figures  6  through  15,  using  the  Colburn-j  modulus,  Fanning  friction 
factor,  and  Reynolds  number. 

Figure  16  Is  a  summary  plot  comparing  the  heat  transfer  and  flow 
friction  characteristics  of  the  matrices  tested. 

Figures  17  and  18  present  "figure  of  merit"  type  curves.   The  j/f 
versus  Reynolds  number  thus  plotted  gives  an  indication  of  the  matrix 
flow  frontal  area  required  for  a  given  pressure  drop. 

Figures  19  and  20  present  curves  of  heat  transfer  power  versus  flow 
friction  power  on  a  unit  area  basis,  evaluated  for  fluid  properties  at 
standard  conditions  of  dry  air  at  500. F  and  one  atmosphere. 

CONCLUSIONS 


Of  the  five  matrices  tested,  the  Solar  No.  2  matrix  with  sinusoidal 
fin  geometry  of  perforated  nickel,  exhibited  the  best  overall  performance 
particularly  at  the  higher  Reynolds  Numbers  tested. 

The  three  stainless  steel  matrices,  Solar  Nos.  3,  A,  and  5,  exhibited 
essentially  the  same  overall  characteristics.   The  triangular  fin  geometry 
of  Solar  No.  5  exhibited  a  slightly  higher  friction  factor  with  a  re- 
sulting lower  1/f  than  the  sinusoidal  fin  geometries. 

The  employment  of  the  cvclic  testing  technique  appears  to  be  an 
effective  means  for  extending  the  range  of  testing  for  the  facility 
originally  designed  for  the  single-blow  transient  testing  technique. 


12 


REFERENCES 

1.  Bannon,  J.  M. ,  Piersall,  C.  H.  Jr.,  Puccl ,  P.  F. ,  "Heat  Transfer 
and  Flow  Friction  Characteristics  of  Perforated  Nickel  Plate-Fin 
Type  Heat  Transfer  Surfaces."  Technical  Report  No.  52,  U.  S. 
Naval  Postgraduate  School,  Monterey,  California,  1965. 

2.  Pucci,  P.  F. ,  Howard,  C.  P.,  and  Piersall,  C.  H.  Jr.,  "The  Single- 
Blow  Transient  Testing  Technique  for  Compact  Heat  Exchanger  Surfaces", 
Trans.  ASME,  Journal  of  Engineering  for  Power,  Vol.  89,  Ser.  A,  No.  1 
Jan.  1967,  pp.  29-40. 

3.  Ball,  Stuart  F. ,  "Experimental  Determination  of  Heat  Transfer  and 
Flow  Friction  Characteristics  for  Several  Plate  Fin  Type  Compact 
Heat  Exchanger  Surfaces."  Master's  Degree  Thesis,  U.  S.  Naval 
Postgraduate  School,  Monterey,  California,  1966. 

4.  Howard,  C.  P.,  "The  Single  Blow  Problem  Including  the  Effects  of 
Longitudinal  Conduction",  ASME  paper  number  64-GTP-ll,  1964. 

5.  Bell,  J.  C.  and  Katz ,  E.  F. ,  "A  Method  for  Measuring  Surface  Heat 
Transfer  Using  Cyclic  Temperature  Variations."  Heat  Transfer  and 
Fluid  Mechanics  Institute,  Berkeley,  California,  1949. 

6.  Locke,  G.  L.  "Heat  Transfer  and  Flow  Friction  Characteristics  of 
Porous  Solids",  TR.  No.  10,  Department  of  Mechanical  Engineering, 
Stanford  University,  Stanford,  California,  1  June  1950. 

7.  American  Society  of  Mechanical  Engineers.   Power  Test  Code  Supple- 
ments, Instruments  and  Apparatus,  PTC- 19. 5 ;4,  Flow  Measurements, 
Chapt.  4,  New  York:   ASME  1959. 

8.  Murdock,  J.  W. ,  "Tables  for  Interpolation  and  Extrapolation  of  ASME 
Coefficients  for  Square-Edged  Concentric  Orifices",  ASME,  paper  No. 
64-WA/FM-6 . 

9.  Kays,  W.  M. ,  "Loss  Coefficients  for  Abrupt  Changes  in  Flow  Cross 
Section  with  Low  Reynolds  Number  Flow  in  Single  and  Multiple-Tube 
Systems,  Transactions,  ASME,  Vol.  72,  pp  1067-1074,  1950. 

10.  Kays,  W.  M.  and  London,  A.  L. ,  Compact  Heat  Exchangers.   Second 
Edition,  New  York:  McGraw-Hill  Book  Company,  Inc.,  1964. 

11.  Davton,  R.  ,T.  ,  et.  al.   "Improved  Measurements  of  Surface  Heat 
Transfer  bv  the  Methods  of  Cyclic  Temperature  Variations,"  BMI-747, 
Battelle  Memorial  Institute,  28  May  1952. 

12.  Chastain,  J.  W.  and  Dinger,  D.A. ,  "A  Study  of  Error  Effects  in 
Measuring  Cvclic-Temperature  Heat-Transfer  Coefficients."   Rpt.  No. 
BMI-1167,  ENGINEERING  (TID-4500,  12th  Ed.),  Battelle  Memorial 
Institute,  Columbus,  Ohio,  15  February  1957. 


13 


w 

0 

a 


■H 

Q 
O 


o 

CO 


-H 


=0 


u 
o 
to 

(0 

<u 

I 

o 

o 
Eh 


14 


e 

Q) 
-H 
(0 

I 

o 


o 


to 

-M 

u 

s. 

4 


0) 

S 

•H 


15 


iWm 


prsnm 


Figure  3.  Heater  and  Test  Section. 


Figure  4.  Test  Section  with  Matrix  Removed. 
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Figure  5.   Nichrome  Wire  Heater. 
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SOLAR  #  1 

Matrix  Material 

Nickel 

Specific  Heat,  cs 

(Btu/lbm  deg  F) 

.106 

Thermal  Conductivity,  k,. 

(Btu/hr  fr  deg  F) 

38.7 

Material  Thickness, 

as 

(in) 

.005 

Total  Heat  Transfer 

Area, 

A    (sq  ft) 

12.768 

Frontal  Area,  A* 

(sq  ft) 

.07022 

Conduction  Area,  As 

(sq  ft) 

.01078 

Free  Flow  Area,  Ac 

(sq  ft) 

.05944 

Matrix  Flow  Length, 

Lm 

(ft) 

.24917 

Matrix  Volume ,  Vm 

(cu  ft) 

.0175 

Matrix  Density,  pm 

(lbm/cu  ft) 

86.1378 

Compactness,  8 

(sq  ft/cu  ft) 

729.7026 

Porosity,  p 

.8464 

Hydraulic  Diameter, 

Dh 

(ft) 

.0046395 

Figure  6.  Geometric  and  Physical  Properties 
of  Solar  #1  Matrix 
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Figure  7.   Heat  Transfer  and  Flow  Friction  Characteristics 
of  Solar  No.  1  Matrix. 
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Matrix  Material 

Specific  Heat,  c 

Thermal  Conductivity,  ks 

Splitter  Thickness,  as 

Total  Heat  Transfer  Area,  A 

Frontal  Area,  Afr 

Conduction  Area,  As 

Free  Flow  Area,  Ac 

Matrix  Flow  Length, 

Matrix  Volume,  V„ 

Matrix  Density, 

Compactness ,  6 

Porosity,  p 

Hydraulic  Diameter,  D^ 

Perforated  Fin  Material 


^ 


"m 
Pm 


SOLAR  if   2 

Nickel 

(Btu/lbm  deg  F) 

.106 

(Btu/hr  ft  deg  F) 

38.7 

(in) 

.005 

(sq  ft) 

12.093 

(sq  ft) 

.07022 

(sq  ft) 

.01034 

(sq  ft) 

.05987 

(ft) 

.25583 

(cu  ft) 

.01797 

(lbm/cu  ft) 

70.3099 

(sq  ft/cu  ft) 

673.143 

.85281 

(ft) 

.0046272 

.0047"  thick 

80/20T 

Figure  8.  Geometric  and  Physical  Properties 
of  Solar  #2  Matrix 
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Figure  9.  Heat  Transfer  and  Flow  Friction  Characteristics 
of  Solar  No.  2  Matrix. 
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SOLAR  #  3 

Matrix  Material 

430 

Stainless 

Steel 

Specific  Heat,  cs 

(Btu/lbm  deg  F) 

.11 

Thermal  Conductivity,  ks 

(Btu/hr  ft 

deg  F) 

12.8 

Material  Thickness,  as 

(In) 

.005 

Total  Heat  Transfer  Area,  A 

(sq  ft) 

14.3588 

Frontal  Area,  Afr 

(sq  ft) 

.06803 

Conduction  Area,  Ag 

(sq  ft) 

.01325 

Free  Flow  Area,  Ac 

(sq  ft) 

.05478 

Matrix  Flow  Length,  1^ 

(ft) 

.24833 

Matrix  Volume,  Vffl 

(cu  ft) 

.01689 

Matrix  Density,  pm 

(lbm/cu  ft] 

I 

94.2656 

Compactness ,  8 

(sq  ft/cu  J 

ft) 

849.8962 

Porosity,  p 

.80517 

Hydraulic  Diameter,  D, 

(ft) 

.0037895 

Figure  10.  Geometric  and  Physical  Properties 
For  Solar  #3  Matrix 
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Figure  11.   Heat  Transfer  and  Flow  Friction  Characteristics 
of  Solar  No.  3  Matrix. 
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SOLAR  #  4 

Matrix  Material 

430  stainless 

steel 

Specific  Heat,  c 

(Btu/lbm  deg  F) 

.11 

Thermal  Conductivity,  kg 

(Btu/hr  ft  deg  F) 

12.8 

Material  Thickness ,  a8 

(in) 

.005 

Total  Heat  Transfer  Area,  A 

(sq  ft) 

14.27694 

Frontal  Area,  Afr 

(sq  ft) 

.06713 

Conduction  Area,  A 

(sq  ft) 

.01229 

Free  Flow  Area,  A<. 

(sq  ft) 

.05484 

Matrix  Flow  Length,  L^ 

(ft) 

.24417 

Matrix  Volume,  Vm 

(cu  ft) 

.01638 

Matrix  Density,  pm 

(lbm/cu  ft) 

88.5918 

Compactness ,  B 

(sq  ft/cu  ft) 

P870.969 

Porosity,  p 

.81690 

Hydraulic  Diameter,  D^ 

(ft) 

.00375168 

Figure  12.  Geometric  and  Physical  Property 
of  Solar  #4  Matrix 
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Figure  13.   Heat  Transfer  and  Flow  Friction  Characteristics 
of  Solar  No.  4  Matrix. 
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Matrix  Material 
Specific  Heat,  cg 
Thermal  Conductivity,  k 
Material  Thickness,  as 
Total  Heat  Transfer  Area,  A 
Frontal  Area,  Ajr 
Conduction  Area,  Ag 
Free  Flow  Area,  Ac 
Matrix  Flow  Length,  1^ 
Matrix  Volume ,  V 
Matrix  Density,  pr 
Compactness,  B 
Porosity,  p 
Hydraulic  Diameter,  D, 


m 


SOLAR  #  5 

430  Stainless 
(Btu/lbm  deg  F) 
(Btu/hr  ft  deg  F) 
(In) 
(sq  ft) 
(sq  ft) 
(sq  ft) 
(sq  ft) 
(ft) 
(cu  ft) 
(lbm/cu  ft) 
(sq  ft/cu  ft) 

(ft) 


steel 
.11 

12.8 
.005 

14.4234 
.06781 
.01148 
.05633 
.26417 
.0179 

81.9345 
805.13313 
.83066 
.0041269 


Figure  14.   Geometric  and  Physical  Properties 
of  Solar  #5 
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Figure  15.   Heat  Transfer  and  Flow  Friction  Characteristics 
of  Solar  No.  5  Matrix. 
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Figure  16.   Comparison  of  Characteristics. 
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Figure  17.   Flow  Area  Goodness  Factors  for 
Solar  Nos.  1  and  2  Matrices. 
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Figure  18.   Flow  Area  Goodness  Factors  for 
Solar  Nos.  3,  4,  and  5  Matrices. 
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Figure  19.  Heat  Transfer  Power  versus  Flow  Friction  Power 
for  Solar  Nos.  1  and  Matrices. 
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Figure  20.  Heat  Transfer  Power  versus  Flow  Friction  Power 
for  Solar  Nos.  3,  4,  and  5  Matrices. 
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Figure  21.      NTU  versus  Amplitude  Ratio. 
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Figure  22.   Perforated  Nickel  Geometry. 
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Figure  23.   Conduction  Parameter  Analog. 
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TABLE  II.      Summary  of  Friction  and  Heat  Transfer  Results 
Solar  #1  Matrix 


FRICTION   RESULTS 


FRICTION 
FACTOR 


REYNOLD'S 
NUMBER 

Nr, 


HEAT  TRANSFER  RESULTS 

Ns  Np   2/3      REYNOLD'S 
NUMBER 


1 


Nr 


NTU 


TESTING 
METHOD 


58553 

35.66 

30134 

69.91 

20595 

106.00 

15680 

136.80 

12574 

170.53 

08867 

254.16 

08029 

267.73 

05971 

394.40 

04418 

523.69 

03663 

642.29 

03713 

660.90 

02736 

912.84 

02419 

1060.99 

02248 

1164.10 

02064 

1293.89 

01995 

1346.47 

01873 

1432.01 

01773 

1557.13 

01652 

1676.73 

01681 

1678.22 

01585 

1854.19 

01486 

2024.35 

01345 

2349.97 

01298 

2504.63 

01265 

2660.83 

01238 

2805.05 

01214 

2945.42 

01202 

3040.39 

01186 

3179.90 

09482 

35.14 

25.57 

S-B 

08829 

68.89 

23.81 

S-B 

05464 

104.47 

14.74 

S-B 

04155 

134.83 

11.21 

S-B 

03475 

168.08 

9.38 

S-B 

02215 

250.50 

5.977 

S-B 

02033 

263.87 

5.48 

S-B 

01401 

388.72 

3.78 

S-B 

00959 

516.15 

2.59 

S-B 

00775 

633.05 

2.09 

S-B 

00812 

660.90 

2.197 

CYC 

00525 

899.71 

1.42 

S-B 

00536 

1060.99 

1.449 

CYC 

00394 

1147.36 

1.06 

S-B 

00400 

1275.26 

1.08 

S-B 

00430 

1346.47 

1.162 

CYC 

00380 

1411.41 

1.03 

S-B 

00378 

1534.71 

1.02 

S-B 

00358 

1676.73 

.966 

CYC 

00326 

1854.19 

.882 

CYC 

00309 

2024.35 

.778 

CYC 

00279 

2349.97 

.754 

CYC 

00283 

2505.63 

.765 

CYC 

00293 

2660.83 

.792 

CYC 

00299 

2805.05 

.807 

CYC 

00305 

2945.42 

.823 

CYC 

00308 

3040.39 

.833 

CYC 

00310 

3179.90 

.837 

CYC 
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TABLE  III.   Summary  of  Friction  and  Heat  Transfer  Results 
Solar  #2  Matrix 


FRICTION  RESULTS 


HEAT  TRANSFER  RESULTS 


FRICTION 
FACTOR 


.53389 
.26494 
.18811 
.15304 
.14937 
.10594 
.08185 
.06003 
.05593 
.04840 
.04383 
.04258 
.03559 
.03521 
.03072 
.02798 
.02720 
.02510 
.03074 
.02505 
.02142 
.02056 
.02078 
.02057 
.02005 
.02021 
.02071 
.02057 
.01848 
.02046 
.02032 
.02019 
.01995 
.01979 
.01951 
.01934 
.01913 
.01890 


REYNOLD'S 
NUMBER 


38.67 

75.32 

110.81 

139.00 

145.58 

212.84 

276.05 

393.56 

421.28 

497.17 

562.64 

570.53 

701.42 

705.01 

823.83 

923.82 

975.36 

1064.41 

1090.38 

1096.74 

1368.82 

1425.27 

1479.07 

1538.31 

1639.80 

1778.11 

1887.96 

2098.59 

2317.65 

2318.17 

2472.89 

2612.20 

2821.31 

2072.19 

3197.52 

3355.29 

3495.91 

3696.46 


2/3 
NStNPr 

REYNOLD'S 
NUMBER 

J 

NR 

.10087 

38.12 

.08623 

74.24 

.05455 

109.21 

.04330 

136.98 

.04233 

143.48 

.02602 

209.78 

.02066 

272.04 

.01317 

393.56 

.01217 

415.24 

.00926 

490.02 

.00966 

562.64 

.00679 

562.26 

.00801 

701.42 

.00598 

694.80 

.00513 

812.02 

.00641 

923.82 

.00428 

961.25 

.00572 

1064.41 

.00531 

1090.38 

.00412 

1096.74 

.00535 

1425.27 

.00545 

1479.07 

.00602 

1639.80 

.00634 

1778.11 

.00658 

1887.96 

.00671 

2098.59 

.00679 

2318.17 

.00688 

2472.89 

.00693 

2612.20 

.00693 

2821.31 

.00690 

2972.19 

.00682 

3197.52 

.00679 

3355.29 

.00674 

3495.91 

.00660 

3696.46 

NTU 


TESTING 
METHOD 


25.60 

S-B 

21.88 

S-B 

13.84 

S-B 

10.98 

S-B 

10.737 

S-B 

6.60 

S-B 

5.24 

S-B 

3.354 

CYC 

3.09 

S-B 

2.35 

S-B 

2.457 

CYC 

1.72 

S-B 

2.038 

CYC 

1.52 

S-B 

1.30 

S-B 

1.630 

CYC 

1.09 

S-B 

1.453 

CYC 

1.349 

CYC 

1.05 

S-B 

1.360 

CYC 

1.385 

CYC 

1.531 

CYC 

1.612 

CYC 

1.672 

CYC 

1.704 

CYC 

1.726 

CYC 

1.748 

CYC 

1.759 

CYC 

1.760 

CYC 

1.751 

CYC 

1.732 

CYC 

1.724 

CYC 

1.710 

CYC 

1.675 

CYC 
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TABLE  IV.   Summary  of  Friction  and  Heat  Transfer  Results 
Solar  #3  Matrix 


FRICTION  RESULTS 


FRICTION 
FACTOR 


REYNOLD ' S 
NUMBER 

Nn 


HEAT  TRANSFER  RESULTS 

Ns  Np  2'3  REYNOLD'S 
NUMBER 


J 


N, 


NTU 


TESTING 
METHOD 


20320 

60.19 

13761 

91.18 

10327 

121.78 

07288 

172.54 

06223 

221.36 

04422 

313.54 

04132 

342.40 

04015 

357.16 

03507 

401.24 

03125 

456.46 

02705 

531.35 

02574 

567.42 

02427 

601.41 

02200 

672.09 

01938 

809.42 

01737 

902.44 

01421 

1144.49 

01378 

1194.66 

01306 

1266.97 

01229 

1380.58 

01193 

1436.51 

01157 

1518.23 

01108 

1675.60 

01123 

1858.27 

00975 

1921.60 

01121 

2013.75 

01123 

2141.30 

01129 

2220.69 

01134 

2320.36 

01134 

2487.33 

01134 

2544.96 

01136 

2658.22 

01132 

2807.06 

01099 

3091.52 

07076 

59.32 

23.30 

S-B 

04325 

89.87 

14.24 

S-B 

03031 

120.03 

9.98 

S-B 

01995 

170.06 

6.57 

S-B 

01464 

218.18 

4.82 

S-B 

01058 

313.54 

3.490 

CYC 

00853 

337.48 

2.81 

S-B 

00930 

357.16 

3.065 

CYC 

00842 

401.24 

2.775 

CYC 

00502 

449.90 

1.65 

S-B 

00634 

531.35 

2.091 

CYC 

00438 

559.26 

1.44 

S-B 

00577 

601.41 

1.901 

CYC 

00379 

662.44 

1.25 

S-B 

00439 

809.42 

1.447 

CYC 

00399 

917.52 

1.316 

CYC 

00346 

1144.49 

1.140 

CYC 

00331 

1266.97 

1.092 

CYC 

00306 

1380.58 

1.009 

CYC 

00285 

1518.23 

.940 

CYC 

00268 

1675.60 

.883 

CYC 

00263 

1858.27 

.866 

CYC 

00269 

2013.75 

.886 

CYC 

00285 

2141.30 

.941 

CYC 

00294 

2220.69 

.970 

CYC 

00306 

2320.36 

1.011 

CYC 

00327 

2487.33 

1.080 

CYC 

00330 

2544.96 

1.089 

CYC 

00331 

2658.22 

1.092 

CYC 

00343 

2807.06 

1.130 

CYC 

00342 

3091.52 

1.126 

CYC 
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TABLE  V.   Summary  of  Friction  and  Heat  Transfer  Results 
Solar  #4  Matrix 


FRICTION  RESULTS 


HEAT  TRANSFER  RESULTS 


FRICTION 

REYNOLD'S 

2/3 
NstNPr 

REYNOLD'S 

NTU 

TESTING 

FACTOR 

NUMBER 

NUMBER 

METHOD 

f 

NR 

J 

Nr 

.37693 

30.98 

.00611 

30.53 

2.00 

S-B 

.20381 

59.13 

.08412 

58.28 

27.52 

S-B 

.14272 

88.76 

.04310 

87.49 

14.10 

S-B 

.11402 

112.87 

.03061 

111.25 

10.01 

S-B 

.07412 

171.76 

.01845 

169.30 

6.04 

S-B 

.05929 

215.99 

.01499 

212.89 

4.74 

S-B 

.06105 

226.98 

.01370 

223.71 

4.481 

S-B 

.03866 

304.36 

.01139 

306.36 

3.732 

CYC 

.04174 

331.25 

.00847 

326.48 

2.772 

S-B 

.03983 

338.77 

.00806 

333.91 

2.64 

S-B 

.03290 

383.18 

.00939 

383.18 

3.078 

CYC 

.03100 

447.77 

.00467 

441.35 

1.53 

S-B 

.03139 

454.43 

.00468 

447.88 

1.532 

S-B 

.02881 

452.28 

.00826 

452.28 

2.708 

CYC 

.02566 

530.33 

.00735 

530.33 

2.409 

CYC 

.02546 

557.45 

.00429 

549.45 

1.40 

S-B 

.02616 

561.71 

.00700 

561.71 

2.292 

CYC 

.02309 

604.56 

.00671 

604.56 

2.200 

CYC 

.02172 

664.38 

.00346 

654.85 

1.13 

S-B 

.02142 

665.20 

.00625 

665.20 

2.048 

CYC 

.01958 

748.02 

.00577 

748.02 

1.890 

CYC 

.01820 

841.71 

.00535 

841.71 

1.751 

CYC 

.01741 

888.19 

.01756 

891.99 

.00511 

891.99 

1.673 

CYC 

.01441 

1120.75 

.00398 

1120.75 

1.303 

CYC 

.01426 

1162.42 

.01291 

1303.44 

.00358 

1303.44 

1.173 

CYC 

.01239 

1415.34 

.01156 

1526.00 

.00324 

1526.00 

1.062 

CYC 

.01103 

1668.92 

.00312 

1668.92 

1.020 

CYC 

.01126 

1807.14 

.00310 

1807.14 

1.015 

CYC 

.01116 

1904.64 

.00313 

1904.64 

1.025 

CYC 

.01161 

1908.14 

.01107 

2041.42 

.00321 

2041.42 

1.072 

CYC 

.01106 

2170.84 

.00327 

2170.84 

1.102 

CYC 

.01156 

2557.98 

.00350 

2557.98 

1.144 

CYC 

.01095 

2711.58 

.00351 

2711.58 

1.148 

CYC 

.01088 

2828.75 

.00356 

2828.75 

1.165 

CYC 

.01086 

2887.46 

.00358 

2887.46 

1.172 

CYC 
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TABLE  VI.      Summary  of  Friction  and  Heat  Transfer  Results 
Solar  #5  Matrix 


FRICTION  RESULTS 


FRICTION 
FACTOR 


REYNOLD'S 
NUMBER 

Nn 


HEAT  TRANSFER  RESULTS 

NStNp  2/3   REYNOLD'S 

NUMBER 


.1 


N, 


NTU 


TESTING 
METHOD 


.22228 

63.05 

.14597 

96.54 

.11431 

121.54 

.08385 

183.46 

.06594 

231.11 

.04962 

337.08 

.04488 

360.72 

.03749 

458.12 

.03439 

482.55 

.03123 

556.72 

.02846 

600.21 

.02472 

708.74 

.02298 

788.74 

.01993 

949.26 

.01907 

964.34 

.01637 

1203.37 

.01548 

1247.65 

.01466 

1405.06 

.01341 

1548.85 

.01320 

1635.14 

.01304 

1804.55 

.01272 

1938.77 

.01207 

2040.53 

.01252 

2061.91 

.01235 

2214.46 

.01233 

2323.94 

.01222 

2480.71 

.01218 

2622.09 

.01212 

2736.70 

.01208 

2876.26 

.01204 

2966.69 

.01198 

3028.94 

.01198 

3069.59 

.06573 
.03924 
.02967 
.01805 
.01393 
.00921 
.00800 
.00670 
.00460 
.00568 
.00426 
.00396 
.00418 
.00355 

.00301 

.00268 

.00250 
.00243 
.00247 

.00249 
.00272 
.00289 
.00319 
.00337 
.00343 
.00353 
.00356 
.00359 
.00361 


62.15 

21.15 

S-B 

95.15 

12.62 

S-B 

119.79 

9.54 

S-B 

180.83 

5.81 

S-B 

227.80 

4.48 

S-B 

337.08 

2.972 

CYC 

355.54 

2.57 

S-B 

458.12 

2.160 

CYC 

475.63 

1.48 

S-B 

556.72 

1.831 

CYC 

591.59 

1.37 

S-B 

698.57 

1.28 

S-B 

788.74 

1.347 

CYC 

946.26 

1.145 

CYC 

1203.37 

.970 

CYC 

1405.06 

.862 

CYC 

1635.14 

.804 

CYC 

1804.55 

.783 

CYC 

1938.77 

.794 

CYC 

2061.91 

.802 

CYC 

2214.46 

.874 

CYC 

2323.94 

.932 

CYC 

2480.71 

1.028 

CYC 

2622.09 

1.086 

CYC 

2736.70 

1.104 

CYC 

2867.26 

1.135 

CYC 

2966.69 

1.145 

CYC 

3028.94 

1.155 

CYC 

3069.59 

1.162 

CYC 
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APPENDIX  A 

The   cyclic  technique  employed  Is   that   described  by  Bell   and 

Katz    (5).      The   assumptions  made   in   the   analysis  are: 

(1)  All   cross   sections   of  the   solid  normal  to  some   fixed  direction 
are   uniform. 

(2)  The  gas  velocity   is   uniform  throughout  the  passage. 

(3)  The  physical  properties   of   the  gas   and  of  the   specimen  are 
constant. 

(A)   The  surface  heat-transfer  coefficient  is   constant. 

(5)    The   inlet   gas   temperature  varies   sinusoidally.      Bell   and 
Katz  have   found  solutions   to   the  problem  under  the   following  three 
alternative   assumptions: 

(a)  The   thermal   conductivity  in   the   solid  is   infinite. 

(b)  The   thermal   conductivity   in   the   solid  is   zero  in   the 
longitudinal    (flow)    direction   and   infinite   in   the   transverse   direction. 

(c)  The   thermal   conductivity  in   the   solid  is   zero  in   the 
longitudinal  direction   and   finite   in   the   transverse   direction. 

Of  the   three   solutions,   that  of   case    (b)    is  most  nearly   satisfied 
in   the   case   of  high  mass    flow  rates  with   its   low  values   of  conduction 
parameter.      A  summary  of  the   theory   for  this   case   is  presented  below: 
Case    (b)   -   Solid  having  infinite  transverse   and   zero   longitudinal 

conductivity. 
Energy  balance: 

Fluid:     mcff^  dx  +  ||£  do)  =  hb(tf  -   ts)    dx  (A-l) 

Solid:      PSASLS  ||^  dx  =  hb(tf  -  ts)    dx  (A-2) 

Inlet   and  Initial  Conditions: 
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Inlet:    tf(O,0)  -  tm  +  A  sin  ojo0  (A-3) 

Initial:   ts(0)  -  tso  (A- 4) 

Initial:  tf(x,0)  *  t0(x)  (A-5) 

where  A  is  the  peak  amplitude  of  the  sinusoidally  varying  tempera- 
ture and  tm  is  some  suitable  mean  temperature.   Introducing  the 
following  dimensionless  parameters: 

(1)  Time       t  »  -^ 

(2)  Position   K   -  f 

to  L 

(3)  Frequency  oj  =  ° 

m  c„ 


(4)  Flow:      M 

WsCsa)      ,  , 
nA 

(5)  Heat  transfer  units  N*u  =  t 


icp 


where : 


u  -   flow  velocity    (ft/sec) 

0  =   time  (sec) 

L  ■  matrix  flow  length  (ft) 

x  =  distance  from  matrix  inlet  in  direction  of  flow  (ft) 

u)0=  angular  frequency  of  temperature  oscillation  (rad/sec). 
Introducing  the  above  dimensionless  parameters  into  (A-l)  through 
(A-5)  and  rearranging  yields: 

Fluid:  f^£   ^£    Ntu(tf  -  t.)  -  0 

6£  ^  6t 

Solid:   6_tg_   M  u  Ntu(tf  -  ts)  =  0 

6t 

Inlet:   tf(0,Y)  »  tso  +  A  sin  cot 
Initial:  ts(0)  =  tso 
Initial:  tf(£,0)  =  t0(0 
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This   system  has   the  periodic  solution: 

tf(F,,T)   =   tm  +  Ae"a?  sin(a)T  +  flO 

ts(C,x)    -   tm  +  Aae"a?[M  Ntu  sin    (ojt   -  QK)    - 

COS    (tOT   -   ft£)l 

where : 

a  =  Ntu/fl  +   (M  Ntu)2] 

a  =  M(Ntu)2/[l  +   (M  Ntu)2]    +  u 

The   ratio  of  the   inlet   temperature   amplitude   to   the  exit   temperature 

amplitude   is 

r       Ntu J 

R  =   e     [  1  +   (M  Ntu)zJ  (A_6) 

This  relationshin  is  plotted  in  Figure  25  as  Ntu  versus  R  for  various 

values  of  M. 

The  solutions  for  the  ratio  of  temperature  amplitudes  for  cases 

(a)  and  (c)  are  presented  below: 

Case  (a).   Infinite  conductivity  in  solid. 

I 


f  -2NtUj_  M2,.    -Ntu,2l2 

r  m     e +  M  (1  -  e    ) 

R   L   1  +  M2   1  -  e-Ntu)^J 


This  solution  is  an  approximation  to  a  Taylor  series  around  to  =  0  when 
NTU  <  2.0. 

Case  (c) .   Solid  having  finite  transverse  and  zero  longitudinal 

conductivity. 

D  _f        Ntu J".  o2M  Ntu    (M2Ntu2-l)  1 V 

R  =  expV  1+(M  Htu)2     L  l  +       6(1  +  M2  Ntu2)  +  "'     J/ 

where : 

o2  ■  w0as  Cs  Ds/2ky 

as  =  plate  thickness 

ky  *   thermal    conductivity  in  transverse   direction 
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This  solution  coincides  with  the  value  found  for  R  for  case  (b) 
whenever  the  series  In  o  can  be  approximated  satisfactorily  by  its  first 
term.   For  a  more  detailed  explanation  on  the  solutions  for  the  preceding 
three  cases  the  reader  is  referred  to  reference  (11).   For  the  treatment 
of  the  three  cases  utilizing  a  liquid  as  the  fluid  the  paper  written  by 
Chastain  (12)  should  be  consulted. 

On  the  use  of  the  three  solutions,  Bell  and  Katz  have  shown  that 
the  ratio  of  the  two  amplitude  ratios  found  for  cases  (a)  and  (b)  are 
coincident  when  M  <  0.1  and  NTU  «  1.0,   In  the  experimental  work  per- 
formed by  Bell  and  Katz,  the  ratio  of  the  two  amplitude  ratios  determined 
for  cases  (a)  and  (c)  was  found  to  be  within  1%  of  unity.   Thus,  if  the 
experimental  values  of  M  are  selected  near  0.1,  and  the  resulting  NTU's 
near  unity,  no  serious  error  is  introduced  by  assuming  zero  thermal 
conductivity  in  the  longitudinal  (flow)  direction,  and  employing  the 
solution  for  case  (b). 
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APPENDIX  B 


Data  Reduction  Relationships 

The  relationships  presented  in  this  appendix  were  used  to  calculate 
the  heat  transfer  and  flow  friction  characteristics  of  the  matrices 
tested. 
Geometry 

The  three  most  significant  geometric  parameters  that  provide  a 
means  of  comparison  between  matrices  are  so  defined  that  there  exists  a 
unique  relation  between  them.   Determination  of  any  two  of  the  following 
fixes  the  third: 

(1)  Hydraulic  diameter 


_    ,     4  x  free  flow  area    .Ac  . 

j),    ss   44  ru   =  —————— ——^——  =   i\        l 

»  h       heat   transfer  area         A 


(B-l) 


(2)  Porosity 


P  = 


free  flow  area    Ac 


frontal  area 
(3)  Area  Compactness 


(B-2) 


8  = 


total 
transfer  surface  area 


AfrL 


total  volume 
If  (B-2)  is  divided  by  (B-3)  the  previous  interdependency  results 

_  P_ 


(B-3) 


rh  * 


B 


(B-4) 

Maximum  Slope 

The  maximum  slope  of  the  generalized  cooling  curve  is  a  function 
of  NTU  and  X  [41. 


Vtfi  -  Hj 


-  ^  (Ntu,  X) 


ma.* 


where : 
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x  =  generalized  time  variable,  tt  =  initial  matrix  temp. 
t  =  hA 


WsCs 


0 


Ntu  =  -. — 

mcr 


tfi  =   air  temp,  upstream  of 

matrix,   0  >  0 

tf2  =   a^r  temp,  downstream  of 
matrix 


Therefore : 


mcr 


Ntu 


TTs^s 


Furthermore 


(B-5) 


tfj  -  t 


i  _  tfo  -  tfi 


tfi  "  M    tfi  -  t± 


+  1 


and 


'(H^-fehns)'  ^-^ 


(E-6) 


combining  (r>-5)  and  (B-6) 


V  tfi  -  tjj  =  v 

d  /  t  \      m 
[  Ntu  / 


_J #  d(tf2-tfi) 

tfi  "  t£       d0 


(n-7) 


considering  the  following  sketch  of  a  generalized  cooling  curve 


0  ,  Time  ■+> 


AA 


d(tf2   -   tf)  I  Y 

dG  IX 


=      dG 


chart  speed 
d(t£2   "   tfi)   =     Y 

(ti  -  tfl)   =   Z 
combining  the  above  measured  and  recorded  quantities  with 

WsCs   matrix  heat  capacity   Cs    -1 

—         s:  ■  ■    55  — —  sec 

mc  fluid  heat   capacity  Cf 

and   equation    (B-7)  ,   yields: 


(! 


f2  -  t; 


fl  "  ti 


(Ntu  ) 


=  Cs   1  Y 
Cf  7,     X 


(3-8) 


l/e  Ac 

with  this  value  of  maximum  slope  and  X  =  -; 

mcp 

enter  Table  1  or  Figure  3  of  Ref.  [21  to  obtain  the  corresponding 

value  of  NTU. 

Cyclic  Relationships 

The  cyclic  tests  were  made  at  the  higher  mass  flow  rates  where 
case  (b)  of  Appendix  A  is  most  nearly  satisfied.   For  this  case 


R 


r   Ntu I 

=  eXpL"  1  +  (M  Ntu)^J 


(B-9) 


Solving  (B-9)  for  NTU,  we  obtain  two  possible  solutions 

=  1  M§  -  (2  M  In  1/R)~2" 
2  M2  In  1/R 


(B-10) 
As  discussed  by  Bell  and  katz  [5] ,  only  the  smaller  magnitude  solution 


has  physical  meaning.   Hence 

1   -  Vl   -    (2M  In  1/r")7 


NTU  = 


(B-ll) 


2  M2  In  1/R 

A  plot  of  equation  (B-10)  is  shown  in  Figure  25.   The  solid  line  portion 
of  the  constant  M  curves  represent  equation  (B-ll). 


A5 


Mass  Flow  Rate 

The  mass  flow  rate,  m,  was  calculated  in  accordance  with  refer- 
ences (7)  and  (**)  : 


where 


m  =   359Kd^FaYVAP  y        (Ibm/hr)  (B-12) 

o  o 


K  =   flow   coefficient,   including  velocity  of  approach    factor 

_C 

=/j   T^    with  C  =  orifice  coefficient  of  discharge  (20) 

and  ~  =  ratio  of  orifice  diameter  to  internal  pipe  diameter, 

d  =  orifice  diameter,  inches 

Fa  =  thermal  expansion  factor  for  orifice  plate 

V  =  expansion  factor  of  fluid 

Y  =  specific  weight  of  fluid  flowing  =  — rzr 

? 

with   P  =   absolute   static  pressure   at   orifice    (lb/ft") 


R  =   53.3    (ft-lhf/lh   -*R) ,   universal  gas   constant    for  air 

L        m 

T  =    [t0C0F)   +  £59.7]    °R      ,    abs.    temp   at  orifice 

^     =  oressure    drop    across    the   orifice    (in.   11-0) 

C  p 

making   the   substitution   of  K  =  ==z~       and  y  =  — rrr  ,    since   the 

fT^F  T 

magnitude  of  the  local  acceleration  of  gravity,  g,  is  taken  as  being 
equal  to  the  magnitude  of  standard  acceleration,  gc,  into  (B-12)  yields: 

**3SV?Tp     dc^Tw(FaY)  cb-13) 

From  reference  (7) ,  Figure  40A 

V  =  1  -  (0.41  +  0.35  Bk)   f 

k 

where  k  =   ratio   of   specific  heats    for  air  ■   1.4 

m  A?n(ln.   Hg   abs.) 
P    (in.   Ilg   abs.) 

and  from  Figure  38  of  same  reference 

Fa  =  1.0 
also:  46 


P°  2 

P  -    (Patm  -  7376   )    .4912  x   144   (lbf/ft  ) 

where : 

Patm  is  the  local  atmospheric  pressure  (in.  Hg) 
Po  is  the  static  pressure  upstream  of  the  orifice  plate  (in.  Hg) 
Inserting  into  (B-13)  the  required  numerical  constants  to  make  the  equa- 
tion dimensional ly  consistent  and  the  above  defined  values  or  relation- 
ships for  P,  T,  Y,  x,  k,  R  and  Fa,  (B-13)  becomes:  Po  -|i 


C       ?  ,   APo   rAPo(Patm  -  13.6 

m  =  589.81^ do  [1  -  (.41  +  .35B4)  pi ]  I  to  +  459.7 

/l  -  e*  Patm  -  Y^-g-  99.02   L 

Reynold's  Number 

By  definition: 

%=  2jS_  (B.15) 

where  G  is  the  mass  flow  velocity  based  on  the  free  flow  area,  Ac 

G  =  "        " 


Ac      p  Afr  (B-16) 

substituting  (B-15)  into  (B-16): 

NR  -  4r-  •  2L  from  (B  -  4)  5l  -  I 


therefore 


K   uAfr    p  p   B 


N*  =(^7T  >  (  6  >  (B-17) 


For  ease  of  calculation,  substitute  (B-3)  ,  B  = 


A 


AfrL 
into  (B-17) 

w   -   ^mL 

No  ■ 


R     Ay  (B-18) 
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Farming   Friction  Factor 

The   following   sketch  depicts   the   flow  system   to  be   considered: 
la  b 

|  ,.     izizz  zzzzzzzzzzzzzi 


2 


&ZZZZZZZZZZZZZZZ2 


va////,  TZZZZZZZ&l 


la  b  2 

For  a  gas    flow  heat  exchanger,   the  pressure   changes    from  section 

1   to  a   and   from  b   to   2   are  very  small   relative   to  the   total  pressure, 
therefore   the   specific  volumes   va~v^,    and  v^~V2.      Since   testing   is 
performed  with   air  at  moderate  pressure,    the   perfect  gas   law  is   applicable. 

The   relation   for  the   flow  stream  pressure   drop   calculation   for  most 
heat   exchanger  cores   is    (10) 


AP 


2gc 


vi         (Kc  +  1  -  p2)  +  2(22.- 

!>!      L  Vl 


l)  +  ff    3E 
Ac     Vi 


entrance 
effect 


flow 
acceleration 

-   (1   -  p2  -  Ke) 


core   friction 


V2  1 
VlJ 


(B-19) 


exit  effect 


Solving    (B-19)    for   f   in   the   isothermal   case   and   recalling   from   (B-l) 

A  L  1 

that  —  =  —       also  v  =  — 

where  pm  =  1  — 2_  with  the  upstream  values  subscripted  by  1  and  the 
downstream  by  2.   Kc  is  the  entrance  coefficient  and  Ke  the  exit 
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coefficient,  both  depend  on  porosity,  geometry  of  the  flow  cross  section, 
and  matrix  Reynolds'  number  [9].   Since  the  ultimate  method  of  calculating 
all  values  for  this  report  was  by  digital  computer,  the  analytical 
expressions  defined  and  formulated  by  Kays  in  reference  [9]  were  utilized. 
(1)  Contraction  ratio,  Cc  p  0.611  +  .045p  +  .344  p  5'7       (B-21) 
(2)  Velocity  distribution  factor,  Kd 
(a)  circular  tubes 
laminar  flow 
Kd  =  1.333 
turbulent  flow 


f  =  0.049  Nr""»2  (Fanning  friction  factor,  circular  tubes) 

Kd  =  1.09068  (4f)  +  0.0588^4?  +  1  (B-22) 

(b)  gap,  laminar  and  turbulent 

Kd  (gap)  -  1  =  0.6 
Kd  (circular)  -  1 

(c)  square,  laminar  and  turbulent 

Kd  (square)  -  1  =  1.17 
Kd  (circular)  -  1 

(d)  triangle,  laminar  and  turbulent 

Kd  (triangle)  -  1  =  1.29 
Kd  (circular)  -  1 

When  utilizing  the  appropriate  Kd  for  the  core  under  evaluation  and  for 

the  existing  type  of  flow, 

Exit  coefficient   Ke  =  1  -  2Kdp  +  p2  (B-23) 

r  <         „.  4      „    1  -  2Cc  +  Cc(2Kd-l)  ,  nix 

Entrance  coefficient:  Kc  ■  9^ (B-24) 

Cc 
Considering  an  order  of  magnitude  approximation,  the  first  term  in  the 
brackets  in  (B-20)  is  the  greatest  contributor  to  the  friction  factor  for 
small  pressure  differentials. 
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Therefore,  the  approximation  that 

Pi   +  P2 

- «  PmaP^Pj  reduces    (B-20)   to 


the   following 

f  -  [2gcpm  §  -    (Kc  +  Ke)    -  £  (Up*)]  |H  <B"25> 

It   is   evident   from   (B-25)    that: 

f«2gc  pm    £7    p-  ,   4rh  =  p/8  and  G  -  m/pAfr 
f*2gc  p,  «(££)(£) 

therefore   for  a  given  matrix 

fa  p3/8 
Colburn  Factor 

Colburn-j    factor  is   defined  as: 

j   =  Nst  NpJk  =  Jl_  Npr%  (B_26) 

V  bCp 

Substituting    (B-16)    for  G  and  multiplying  —yields: 


.        hA       Ac     .,      &        .         „_..  hA 

j   ■ t-  .  NDr<3        but  NTU  -     — 

mCr,     A         ^l      '  mc 


therefore: 

Ac         34 

j   =  NTU  —  Npy7*    combining  equations    (B-l) 

and    (B-A) 


Ac  _  p 
A      !   8L 


thus 


j  =NTUNpr*  1. 1  (B-27) 

It  is  apparent  from  (B-27) 

J    8 
Heat  Transfer  Power  and  Flow  Friction  Power 

An  evaluation  of  the  heat  transfer  power  versus  flow  friction  power 
is  of  interest  in  that  it  is  a  measure  of  relative  performance. 
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Beginning  with   the  definition  of   the  Colburn-j    factor: 


J   -  NstNpr 

Nst-     N|lu       = 
NrNpr 

hDh 
k 

k 

1 

% 

therefore : 

Cpy  NR  pr 


solving   for  h 

N    Vj      Dh  R  J 


P1 


If  Cp    ,   y    ,   and  Npr  are  evaluated  at   standard  conditions   of  dry  air 
at  500°F  and  one   atmosphere    [10] 
cp  =  0.2477  RTU/lb-F 
M  =   0.0678   lb/hr-ft 
p  =  0.0413   lb/ft3 
Npr  =  0.671 
Equation   (B-28)   becomes: 

1  RTIT 

hstd  -  0.02195(^)    QU)  ggffig   _   .F  (B_M) 

The   flow  friction  power  per  unit   area  is   derived  as   follows: 

AP 


f  =  2gc  pm  gy-  /  Bl)  rearranging 


E  =  AP  Bi 
L 

therefore  : 


AP     lb.  =  ~i ,   since  G  -  pmV 

L  2gc  pm 


E=2i^      pffeVV     <"*  (B"30) 


3 
m 

When  equation    (B-30)   is   evaluated  at  standard  conditions    (y   and 

pm  values   given   above) ,   the   flow   friction  power  per  unit   area  at   Standard 

conditions    is: 
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ESTD     -  1.11  x  ID"'       j|-       f    (MJfr    „P/ft*  (b-31) 

For  comparison  purposes   the  surface  geometries  were  reduced  to  a 
common  hydraulic  diameter  of  11  ■  2.   x   10     ft. 
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APPENDIX  C 

Description  of  Equipment 

The  equipment  utilized  to  obtain  the  data  from  which  the  heat 
transfer  and  flow  friction  characteristics  were  determined  was  basically 
the  same  for  both  methods  of  testing.   Only  one  change  was  required,  and 
that  was  in  the  power  equipment  to  produce  the  necessary  temperature 
change  corresponding  to  the  method  employed. 

This  appendix  will  be  divided  into  two  main  subsections; 

(1)  Equipment  common  to  the  cyclic  and  single  blow  testing 
systems 

(2)  Power  generating  equipment  utilized 

1.   Equipment  common  to  the  cyclic  and  single-blow  testing  systems 
Fluid  Source; 

Air,  the  working  fluid,  was  drawn  through  the  test  apparatus  by 
a  30HP,  multistage,  Spencer  Turbo-Compressor  rated  at  550cfm  operating 
on  a  220  volt  a.c.  power  supply. 
Flow  Metering  System: 

Flow  rates  were  measured  by  means  of  an  ASME  standard  orifice  sec- 
tion designed  for  convenient  changing  of  orifice  plates  [7].   Thin  plate 
concentric  orifices  with  throat  diameters  from  .775  to  2.31  inches  were 
employed.   The  inside  pipe  diameter,  d,  was  3.08  inches.   Pressure  taps 
were  located  one  pipe  diameter  upstream  and  one-half  pipe  diameter  down- 
stream from  the  orifice. 
Fluid  Heater  System: 

The  air  heater  section  consisted  of  28  nichrome  .0031"  diameter 
wire  heaters.   Nichrome  was  selected  because  of  its  high  resistivity,  low 
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thermal  conductivity,  and  low  specific  heat. 

The  nichrome  wire  was  wound  on  bakelite  frames,  two  heaters  to  a 
frame.   The  wires  were  wound  1/32"  apart,  50  and  52  wires  respectively 
to  each  heater,  and  were  formed  into  a  rectangular  mesh.   The  heaters 
on  each  frame  were  electrically  connected  in  parallel  and  then  connected 
through  a  toggle  switch  to  voltage  bus  bars. 

A  schematic  wiring  diagram  for  one  frame  is  shown  below: 
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All  heaters  selected  operate  in  parallel  when  power  was  applied.   The 
number  of  heaters  used  was  dependent  on  the  testing  method.   For  single 
blow  testing  the  number  of  heaters  used  was  determined  by  the  number  re- 
quired to  effect  a  20°F  rise  for  a  given  mass  flow  rate  when  a  220  volt 
source  was  applied,  e.g.,  at  m  =  790.  lbm/hr  and  V  ■  220  a.c,  22  heaters 
were  required.   For  cyclic  testing  all  28  heaters  were  used  while  the 
input  voltage  was  varied. 
Matrix  Holder  and  Test  Section  (Figure  4) 

The  matrix  holder  and  test  section  were  constructed  from  polyethylene 
plastic.   The  test  section  casing  contained  a  removable  frame  on  which  a 
set  of  five  thermocouples  connected  in  series  was  mounted,  T2.   It  also 
contained  two  static  pressure  taps,  one  upstream  and  one  downstream  from 
the  matrix  holder.  The  components  making  up  the  casing  and  holder  were 
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machined  to  close  tolerances,  thereby  minimizing  air  leakage,  and  insur- 
ing good  alignment  of  the  flow  channel  through  the  system.   The  flow 
channel  was  3-3/16"  by  3-3/16".   Matrices  with  flow  lengths  up  to  3" 
could  be  tested.   Another  removable  thermocouple  frame  containing  a  set 
of  five  thermocouples  connected  in  series,  T^  ,  was  mounted  in  the  exit 
plane  of  the  matrix  holder.  The  matrices  were  placed  in  the  holder  and 
held  in  place  by  insertion  of  styrofoam  plastic  on  four  sides  of  the 
matrix  to  insure  a  snug  fit  in  the  holder  and  proper  position  relative 
to  the  gas  flow. 
Pressure  Measuring  System 

Pressure  taps  were  located  upstream  and  downstream  from  the  matrix 
holder  and  on  either  side  of  the  orifice.   Each  tap  was  connected  to  an 
appropriate  manometer  or  draft  gage  via  Imperial  Company  "poly-flow" 
tubing.   The  following  instruments  were  utilized: 

(1)  E.  Vernon  Hill  and  Company  Type  "C"  Micromanometer,  0-125" 

(2)  Ellison  Draft  Gage  Company,  0-3.0"  inclined  gage 

(3)  Ellison  Draft  Gage  Company,  10"  manometer 

(4)  Ellison  Draft  Gage  Company,  20"  manometer 

(5)  Merriman  Instrument  Company,  50"  manometer 

(6)  Merriman  Instrument  Company,  120"  manometer 

Any  one  of  the  above  instruments  or  combination  could  be  used  to  measure 
the  differential  pressures  at  the  orifice  or  matrix.   Cross  checking  of 
the  various  instruments  assured  reliable  operation. 
Temperature  Measuring  System:  (Figure  1) 

All  temperatures  were  measured  with  thermocouples.   There  were  four 
locations  in  the  system  where  a  temperature  was  monitored;  prior  to  the 
flow  straightner,  before  and  after  the  matrix,  and  before  the  orifice. 

55 


A  30  gage  copper-cons tantan  thermocouple  upstream  from  the  orifice 
plate  was  referenced  to  an  ice  junction.   The  output  was  read  from  a 
Rubicon  Company  Portable  Precision  Potentimeter. 

Four  sets  of  30  gage  iron-cons tantan  thermocouples  were  utilized 
at  the  remaining  three  positions.   Each  of  the  four  sets  was  made  up 
with  five  thermocouples  connected  in  series.  All  20  thermocouples  were 
of  similar  construction  and  were  fabricated  from  the  same  spools  of 
iron  and  constantan  wire.   Thirty  iron- cons tantan  couples  were  made  up 
and  tested  for  temperature  response  (1)  with  the  aid  of  a  Leeds  and 
Northrop  Company,  Type  K-3  Universal  Potentimeter.   From  the  couples 
tested  the  20  thermocouples  that  had  the  lease  diviation  from  the  average 
output  were  selected.  The  maximum  deviation  from  the  average  output 
of  the  thermocouples  selected  was  018  +  .2  yv. 

Two  of  the  four  sets  of  iron-cons tantan  thermocouples  were  inserted 
into  a  1/4"  diameter  aluminum  tube  which  served  as  a  radiation  shield. 
The  tube  was  mounted  in  a  frame  located  at  the  exit  of  the  contracting 
cone.   All  ten  thermocouples  in  the  two  sets  were  individually  wrapped 
with  teflon  tape  to  prevent  shorting.   Each  set  was  designated  as  Ti. 
Another  set,  T2 ,  of  five  thermocouples  connected  in  series  was  securely 
mounted  in  a  frame  in  the  test  section  located  immediately  upstream  from 
the  matrix.   The  remaining  set  of  thermocouples  was  mounted  in  a  frame 
located  in  the  exit  plane  of  the  matrix  holder.   The  latter  set  was 
designated  T3.   The  output  of  T2  and  T-  were  individually  referenced  to 
the  outputs  of  the  two  T- 's  located  upstream  at  the  entrance  to  the 
system.   Therefore,  the  output  of  (T.  vs  T.)  represented  the  magnitude 
of  the  temperature  rise  of  the  heated  air  above  ambient  entering  the 
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matrix,  and  (T.  vs  TO  sensed  the  temperature  rise  of  the  exiting  gas 
from  the  matrix  above  ambient. 

For  the  single  blow  tests  a  record  of  the  output  of  (T^  vs.  TO 
as  a  function  of  time  was  made  for  each  run.   For  the  cyclic  tests  the 
outputs  of  (T^  vs  T2)  and  (T^  vs  TO  were  recorded.   The  recording 
instrument  was  a  Mosely  Company,  Model  7100B  dual  channel  strip  chart 
recorder.   Sixteen  fixed  full  span  settings  or  any  magnitude  of  voltage 
span  from  1  millivolt  to  100  volts  were  available  for  selection.   Twelve 
fixed  chart  speeds  in  the  range  of  1  inch  per  hour  to  two  inches  per 
second  were  provided. 
2.   Power  Generating  Equipment  Utilized. 

For  the  single  blow  test  runs  the  power  to  the  heaters  was  supplied 
by  a  240  volt  a.c. ,  60  cycle  source  regulated  by  a  General  Radio  Company 
Type  W20HM  "Variac"  auto  transformer,  0-240  volt,  8  amp  load  capacity. 
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With  the  arrangement  shown  above  schematically,  the  amount  of  volt- 
age supplied  to  the  resistive  heaters  could  be  simply  regulated.   The  on- 
off  knife  switch  provided  a  means  by  which  the  power  could  be  quickly 

57 


secured. 

The  cyclic  technique  required  a  means  of  creating  a  low  frequency 
sinusoidal  temperature  wave.   The  frequency  range  was  dictated  by  the 
core  material  and  mass  flow  rate  for  each  test  run.   This  restriction  is 
implied  by  the  requirement  of  keeping  the  flow  parameter,  M  <  0.1.   At 
the  maximum  mass  flow  rate  of  2000.  lbm/hr  used  during  testing  the 
maximum  frequency  required  was  2.1  cycles  per  second. 

The  system  finally  employed  to  fulfill  the  requirements  of  the  cores 
reported  consisted  of  a  1.5  KW,  3  phase,  1725  RPM,  inductively  wound  G.E. 
amplidyne  motor  generator  and  a  Hewlett  Packard,  Model  203A  low  frequency 
function  generator.   The  schematic  drawn  below  indicates  the  manner  in 
which  these  components  were  used. 
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The  function  generator's  voltage,  amplitude,  and  frequency  output 
controlled  the  field  excitation  of  the  amplidyne.   The  output  level  of 
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the  amplidyne  was  raised  above  electrical  ground  by  an  amount  controlled 
by  the  setting  of  a  high  wattage  resistance  box  connected  in  series  with 
a  250  volt  coarsely  regulated  d.c.  generator.  The  magnitude  of  the  d.c. 
voltage  fixed  the  mean  temperature,  tm,  of  the  flowing  gas. 

The  frequency  response  was  relatively  flat  up  to  approximately  1.2 
cps.   Its  peak  to  peak  voltage  output  in  this  frequency  range  was  80 
volts.   As  the  system  was  driven  faster  by  the  function  generator  the 
peak  to  peak  voltage  amplitude  decreased  to  about  30  volts  at  2.5  cps. 
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APPENDIX  D.   PERFORATED  NICKEL  PROPERTIES 

The  perforated  nickel  plate  fins  of  the  Solar  #2  matrix  ia  an  electro- 
deposited  metallic  sheet  of  integral  structure  manufactured  by  Perforated 
Products,  Incorporated,  and  is  designated  by  the  manufacturer  as  80/20T. 
The  material  is  .0047  inches  thick  with  14%  open  area  formed  by  slotted 
holes  nominally  .0023  inches  by  .034  inches,  the  centers  of  which  are 
spaced  .0125  inches  in  the  flow  direction  and  .050  inches  transverse. 

For  a  solid  splitter  plate  or  solid  fin,  the  conduction  parameter  is 

defined  as 

X  =  ks^-s 
L  C 

where  As   is   the   solid  cross-sectional   area  for  conduction  which   is   equal 

to  the  product  of  plate  width   and   thickness,    and  L  is   the   conduction  path 

length  which  is   identical   to  the  matrix   flow  length.      The  grouping  of 

terms  r— -       represents   the  thermal   conduction  resistance,    thus,    the 
^sAs 

conduction  parameter  may  be  expressed  as 

1  t. 

X   =  - — r  where        R 


RthC  th  ksAs 

Using   this    form,   the   conduction  parameter  for  a  perforated  material 
will  be   represented  by 

1 


Xk  = 


Wkc 


It  is  useful  to  define  a  relative  conduction  shape  factor,  S  ,  such 


that 


Rth       *k 


<*th>k  * 

Thus   S     is   seen  to  be   the  ratio  of   the   conduction  parameter  of  a  perforated 

material  to   the  conduction  parameter  of  the  same  material  imperforated. 
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By  analogy  between  conduction  heat  transfer  and  electric  current, 
the  ratio  of  electrical  resistances  to  current  flow  of  two  specific 
geometrical  shapes,  is  equal  to  the  ratio  of  thermal  conduction  resist- 
ances to  heat  flow  in  the  same  direction.   Thus 

r  _  Rth 
rk   (Rth)k 

Two  dimensional  models  were   constructed   from  electrically  con- 
ducting paper    (Teledeltos  paper)    representing   repeatable  sections   of 
the  perforated  material   and  its   unperforated   counterpart.      Figure   26   is 
a  plan  view  of  the  perforated  material,   the   shaded  portion   representing 
the   repeatable   section.      Figure   27  below  it   shows   the   two  shapes    for 
which   the   resistances  were   compared.     Measuring  the   resistances  with 
an   ohm-meter  yielded   the   ratio  of  electrical   resistances 

I_  „  Rth        m   s*  =  Mc 
rk      (Rth)k  "  * 

Finally,   the  conduction  parameter  for  the  perforated  material  is  obtained 

X,  =  S*X       where  X  =  ^^ 
k  L  C 
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